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Abstract 

We present an adequate analytical approach to the description of nonlinear vibra- 
tion with strong energy exchange between weakly coupled oscillators and oscillatory 
chains. The fundamental notion of the limiting phase trajectory (LPT) correspond- 
ing to complete energy exchange is introduced. In certain sense this is an alternative 
to the nonlinear normal mode (NNM) characterized by complete energy conserva- 
tion. Well-known approximations based on NNMs turn out to be valid for the case 
of weak energy exchange, and the proposed approach can be used for the descrip- 
tion of nonlinear processes with strong energy exchange between weakly coupled 
oscillators or oscillatory chains. Such a description is formally similar to that of a 
vibro-impact process and can be considered as starting approximation when dealing 
with other processes with intensive energy transfer. At first we propose a simple 
analytical description of vibrations of nonlinear oscillators. We show that two dy- 
namical transitions occur in the system. First of them corresponds to the bifurcation 
of anti-phase vibrations of oscillators. And the second one is caused by coincidence 
of LPT with separatrix dividing two stable stationary states and leads to qualitative 
change in both phase and temporal behavior of the LPT (in particular, temporal 
dependence of the amplitude becomes resembling that for vibro-impact vibrations). 
Next problem under consideration relates to intensive intermodal exchange in the 
periodic nonlinear systems with finite (n>2) number of degrees of freedom. We con- 
sider two limiting cases. If the number of particles is not large enough, the energy 
exchange between nonlinear normal modes in two-dimensional integral manifolds 
is considered. When the number of the particles increases the energy exchange 
between neighbor integral manifolds becomes important that leads to formation of 
the localized excitations resembling the breathers in the one-dimensional continuum 
media. 
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Introduction 



Understanding the mechanisms of energy transfer and locahzation is one of the key prob- 
lems in physics, mechanics, chemistry and biology. This problem is strongly connected 
with searching the elementary excitations preserving their characteristics under mutual 
interactions. The nonlinear normal mode (NNM) that is an extension of well-known nor- 
mal mode in the linear systems (LNMs) is the important type of elementary excitation. 
Contrary to LNMs, the NNMs do not possess the superposition property. Besides, the 
NNMs can demonst rate a spatial loca l izatio n and their nurnber ma y exceed the number of 
degrees of freedom Manevitch et al.l (1l989l ): IVakakis et al. I (1l996l )]. In the linear systems 
the energy cannot be transferred by LNM (because they have homogeneous structure and 
infinite prolongation). As for some combination of LNMs, in particular, the wave pack- 
ets, they c a n tran sfer the energy, but inevitably decay because of dispersion phenomenon 



Witham. I fll974 )]. If nonlinearity is taken into account, the solitonic mechanism of en- 



ergy transfer (by non-topologic al or topologi cal solitons as well as by envelope solitons 



(breathers)) becomes possible [Scott. I (120031 )]. But the conventional results in this field 
relate to infinite continuum systems or oscillatory chains. In this connection two im- 
portant questions arise: (1) what is a mechanism of energy transfer in relatively small 
periodic oscillatory chains? And (2) how does transition to localized excitations proceeds 
when the number of the particles in periodic chain increases? 

The object of the paper is to give the answers on these questions. At that is shown 
that the conception of Limiting Phase Trajectory (LPT) has a key role in the understand- 
ing the regularities of intensive en ergy transfe r in th e oscillatory chains. This conception 
was first introduced in the paper |Manevitch.l fl2007f )] in the application to two degrees of 
freedom (2DoF) system. We consider first the intensive interaction of degenerate NNMs 
(which may be correlated with q-breathers in the Fourier space of the system) in the 
invariant manifolds, described by LPTs, as mechanism of slow energy transfer by peri- 
odic traveling (envelope) waves. Then it is shown that a closeness of the frequencies, 
corresponding to different integral manifolds, provides formation of stationary or mobile 
localized excitations resembling the breather in the infinite chain when the number of the 
particles increases. 



1 Two weakly coupled nonlinear oscillators 

Let us consider first the simplest nonlinear problem of energy transfer in two weakly 
coupled nonlinear oscillators with cubic restoring forces (FigH]). 
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Figure 1: The system under consideration 



This problem can be described by the following system of two nonlinear equations 
(in dimensionless form): 
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Cl and C3 are the hnear and nonhnear stiffnesses of the first and second oscillators, 
respectively, and C12 is the stiffness of the coupling spring. Introducing the complex 
variables: 
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and slow time T2 = eti (along with the fast time ri), one can use the following 
two-scale expansions 



(n, T2) = ^ ^j,n (n, rs) e", j = 1, 2. 



(4) 



After corresponding calculations taking into account Eqs. ([IM]) we arrive at the 
equations of the principal asymptotic approach 



£ + ^/5/2-3^a|/2|/i = 0, 
^+tf3h-3ta\fi\f2 = 0, 



dT2 



(5) 



whic h describe, as well as the mech anical system discussed, others including optic 
couplers Akhmeriev. Ankiewicz. I (119921 )] . This system (Fig. [1]) is fully integrable and 
has two integrals: 



i^ = /3(/2/r + /i/2)-2«(i/ii' + i/2i') 



iV=l/l|'+l/2|' 



(6) 
(7) 



The best way to address this is to use Eq. [7] and the coordinates and A, where 



/i = ^cosee'^\ /2 = v^cos^e^^^ A = (5i = * 



2, 



4^ = /3 sin A, sin 26^ = 2f3 cos 26 cos A + -aN sin 46. 

dT2 dT2 2 

The integral of Eq. [9] has the form following from Eq. [6] 



(9) 



Hr, = COS A + sin 26, 



where k 



3a N 



a > 



(10) 



The system [9] is strongly nonlinear even in the case of initially linear problem. 
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Figure 2: Phase trajectories in the A) plane for: a) k = 0.2, h) k = 0.4, c) k = 0.55, 
d) k=0.9, e) k = 1., i) k = 1.5 
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Let us present plots of the phase trajectories for different values of k in Fig. [21 
(Because of the phase plane periodicity one only has to consider the two lower quadrants.) 

Let us show that we can consider the LPT as fundamental solution (similarly to 
NNMs) whose behavior determines the second dynamic transition in the behavior of the 
oscillatory system. The LPT, which is far from the stationary points, can then be used 
as a generating solution to construct close trajectories with strong energy transfer. 

The value of H for the LPT in integral [10] is equal to zero. Therefore the variables 
9 and A in this case are connected by the equation 

cos A = -A; sin 2^, (11) 

so that sin A = ±a/i - P sin^ 26'. 

Then the first of the equations [9] can be written as: 



±13^1 - F fAn^2e (12) 



The solution of Eq. f[T2|) for the plus sign is the elliptic Jacobi's function 6 = 
{l/2)am{2j3T2, k). Because > 9 > 7i/2 by definition, one can use the negative sign for 
{2n - l)7r < 2/3t2 < 2mT, n= 1,2,3... 

Thus we arrive at the final solution 

6* = - |am (2/5r2, k) \ , A = ± arccos [k sn {2(3t2, k)] , (13) 

with period K{k), i.e., the complete elliptic integral of the first kind (for the in-phase 
oscillations). 

The solution for the out-of-plane oscillations is 

= - I am (2(]ti, k)\ , A = tt ± arcsin [n (2/3ri, k)] . 

The periodic functions 0131) are not smooth; A(2/5r2) has breaks at the points 2(3t2 = 
(2n — l)7r, n=0,l ,. . . , and 9{2(3t2) has discontinuities in the derivative at these points (in 
terms of distributions d9/dT2 = (2/9/7r)A). Plots of 9{2(3t2) for two values of parameter k 
are presented in Figs. [3]and[H The value k = 0.5 corresponds exactly to the first dynamic 
transition. However, the solution for the LPT (Fig. [3]) is still close to that of the linear 
case, except for a small change of period. Only for values of k that are close to unit the 
deflections from an exact saw tooth profile and the change of period become noticeable. 

The second dynamic transition occurs when k = 1. In this case one can find a simple 
analytical solution corresponding to the LPT: 

/■^^ d (29) 1 1 - 

2Pt2 = / 9 = - arcsin \, , (14) 

^ Jq cos29' 2 1 + ' ^ ' 

It can be seen from Eq. [14] that the LPT actually coincides with a separatrix: if 
k —>■ 1, 9 ^ 7i/A when T2 oo. 

It is convenient to introduce two non-smooth functions t{t2), e(r2) (Fig. [5]). Similar 
functions (but wit h alternating signs of dependent variables) were introduced first in 



Pilipchuk. I (119881 )] to calculate the smooth vibrational regimes close to v ibro-impact 



ones. The development of corresponding techniques is also presented in [Pilipchuk. 



( 119881 )] ■ Here we use similar notation. 
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Figure 3: The function 29 = \am (2x2, k)\ for k=0.5 




Figure 4: The function 26 = lam {2t2, k) \ for k=0.9 




Figure 5: The non-smooth basic functions r(r2/a), e(r2/a), where 2a is a halved period 
(in time T2) 
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We would like to show that the natural area for application of these non-smooth basic 
functions is the description of beats (using the variables 9 and A) and close trajectories 
with strong energy transfer. Actually in the case k = Q (the linearized system) the solution 
(IH]) can be rewritten in the form = (7r/2)r, A = (7r/2)e, r = r(r2/a), e = e(r2/a), where 
a = 71/2(3 (exactly as in a vibro-impact process with velocity A = 7r/2). After introducing 
the basic functions r(r2/a), e(r2/a), we can present the solution as 

e = Xi(r) + Fi(r) e (^) , A = X,ir) + Y,{t) e (^) , (15) 
where the smooth functions Xj(r), Yi^r) satisfy Eqs.([nD: 



Yi 



d_ 

dr 

d_ 

dr ] y 



^a/5[sin(X2 + Y^) T sin(X2 - 1^2)] 



2 



al3[cot2{Xi + Yi) cos(X2 + F2 



^cot2{Xi - Yi) cos(X2 - F2)] + fa X[cos 2(Xi + Yi) t cos 2(Xi - Yi)] 

Then, we can search for the solution of Eq. f|T6|) in the form of power expansions in 
the independent variable r: 

00 00 
X, = Y,X,,it\ Y^ = Y,Y,,it\ 1 = 1,2 (17) 

where the generating solution is the linear beat: Xi o = 0, Xi 1 = 7r/2, Yi q = 0, 
^2,0 = 0, Y2fi = n/2, satisfying exactly the 9 — A equations for the case of the strongest 
beat. It can be proved that the presentation (ITSl) . taking into account (fT71) . actually 
recovers the exact solution of the nonlinear problem for the most intensive energy transfer 
between the oscillators. As this takes place, the expansions (fT7|) restore the exact local 
representation of the corresponding elliptic function (near r = 0), but the expressions 
(IT^ allow the prediction of the exact global behavior of the system. It is important to 
note that, even for large enough values of k, the solution appears close to that of linear 
beats; the only change is the barely seen curvature of the lines that are straight for linear 
beats, and a change of the period. 

One can find corresponding corrections by considering the next order of approxima- 
tions, namely, Xi,o = 0, Xi,i = a(3, Xi,3 = -(2/3)(a/3)3p, 1^2,0 = vr/2, ^2,1 = 7r/2, 
y2,i = 2a(3k, l2,3 = 4/3(a/9)'^, . . . , which coincide with those in the expansions of the 
exact solution. 



2 Resonant energy transfer and energy localization 
in finite FPU systems with periodic boundary con- 
ditions 

FPU systems with periodic boundary conditions 



Let us consider a common FPU-system [Fermi et al. I (Il955l )]. the Hamiltonian of which 
can be written as follows: 
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where qj and pj are the coordinates and conjugate moments, respectively, and N is 
the number of particles. 

The transformation to normal coordinates is a linear canonical transformation 



iV-l 



9i = Yl ^i'k^k 



k=0 



with 



r 1 



Sj,k 
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+ l,...,iV-l 



that allows us to write the quadratic part of Hamilton function as the sum of energies 
of independent oscillators: 



N-l 



H2 = J2W^ + -^l^^ 



(20) 



k=l 



where and 77^ are the amplitudes and momentums of the NNMs. (The coordinate 
corresponding to the center of mass motion (^0) is removed from ( l20l) ). 
The eigenvalues 



ciJfc = 2 sin( 



A; = 0, 



N 



(21) 



are restricted by the top-frequency un/2 = 2. The zero eigenvalue uq = corresponds 
to the motion of mass center. The spectra of periodic FPU chains with various numbers of 
particles are shown in FiglHl The short horizontal lines correspond to the eigen-frequencies 
of the chain with N particles, and all of the eigenvalues in the interval < cj^ < 2 are 
twice degenerate. Let us note that in the case of the system with even number of particles 
the gap between frequencies of the top non-degenerate mode k=N/2 and nearby modes 
decreases fast enough when the number of particles increases. This fact is important for 
further analysis of instability of top-frequency vr-mode. 

One should note that in spite of harmonic part of Hamilton function does not sen- 
sitive to the value of the phase 7 in the canonical transformation (fT9l) . it is not valid 
for the nonlinear interaction. The nonlinear part of the Hamilton function for the pe- 
riodic /3-FPU chai n has the simplest representation with the specific choice of 7 = 7r/4 



Poggi and Ruffo. I (119971 )]: 



A 

SN 



N-l 



^k^l^m^nCk,l,m,n^k^l^m^ni 



(22) 



k,Lm,n=l 



where 
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Ck I m,n ^ 



k+l+m+n ~r '-^k+l—m 



-I-/ — IT). — n 



k—l+m—n ~t~ Afc__;_m_|_n 



with 



^ ^ otherwise. 
The respective equations of motion are written in the form: 



where 



N-l 
Lm,n=l 



(23) 



Figure 6: Spectra of eigenvalues of FPU chain with periodic boundary conditions for 
various number of particles 



If we assume that parameter of non-linearity P is equal to zero, the equations ( ]23l) 
describe the set of independent oscillators. What happens when the parameter j3 is 
not equal to zero? Wh ether some NNMs preserve their independence? It is known 



Poggi and Ruffo. I (119971 )] that the modes with indexes 



k=0, N/4, N/3, N/2, 2N/3, 3N/4, 
are decoupled from other modes and they are one-mode particular solutions of 
eqs.f l23l) . The mode with zero index is trivial. The top-frequency mode with k = N/2, 
which is frequently named as "vr-mode" or "zone-boundary mode", is one of the exact 
solution of equations (1231) . The remaining modes, which are equidistant from the N/2, 
form a series of two-mode invariant manifolds in the Fourier space of the system. It is 
necessary to note that this result is valid for the /3-FPU lattice. 

But what we can say about the other types of elementary excitations? It is clear 
that the linear properties of normal modes have to be approximately preserved if the 
parameter of nonlinearity (3 is small (or, that is the same, the amplitudes of modes are 
small enough). The traditional conception is that the linear oscillators far from resonant 
conditions do not interact. Moreover, the presence of 1:1 resonant conditions (between 
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modes with indexes k and (N-k)) in the p eriodic weakly nonhnear FPU-chain does n ot 



disrupt the integrabihty of original system Rink. I (120011 ) jHenrici and Kappeler. I (120081 )]. 



To analyze a dynamics of the system we introduce the complex amplitude of normal 
mode \l/ as follows: 

in the terms of which the equations ( !23ll one can be written as the coupled nonlinear 
equations: 



i,m,n=l 

We need to extract the mono-frequency motion to concentrate our attention on non- 
linear evolution. Because we consider the processes in the time scale that is essentially 
large than 1/uJk, we have to introduce the time hierarchy: 

To = t, Ti= ETo, T2 = ^Vq, (25) 

where the "fast" time tq corresponds to the natural time scale of the modes and the 
"slow" times Ti and T2 allow to study the behavior of envelope functions. Keeping in mind 
the dependence of chain properties from the chain length, one should consider the value 
1/A^ as a small parameter e. 

Let us introduce the envelope function ipk as follows: 

V^fc = v^(Xfc,i + £Xk,2 + £'^Xk,3 + •••)• (26) 



The multi-scale expansion based on the equations (1251 [26]) allows us to construct the 
dynamical equations for the main-order approximation amplitudes Xk,i- One can show 
that for the /3-FPU system equations ( l24l) admit the particular solutions in the form of 
a pair of degenerate modes with indexes (k, N-k) when the eigenvalues Uk are different 
enough. We call these degenerate modes as "conjugate" ones. 

Interaction of conjugate modes. 

As it was mentioned above, in the small FPU chain the modes, which are not conjugate, 
have the frequencies which are differed significantly. In such a case we consider first the 
interaction between conjugate modes. One can see that the amplitudes of main-order 
approximation Xk,i are independent on the "fast time" Tq. Really, in the lower-order by 
small parameter e we get: 

The equations of the next order by small parameter e are read as follows: 
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The appropriate solutions of equations flTTI) are 



(9to 9ri 



XAr-fc,2 = 0, Q^^'^ = 



So the equations of the next order by small parameter are: 



5/2 . KsXN-kXl — Rk /^r 
^ ^^^^ J + l/^llXAf-A:,!! +'t2|Xfc,l| )XN~k,l + 



The right hand side items Rj in equations (!28|) contain both the "resonant" and 
the "non-resonant" parts of the interactions that dependent on the fast time tq. The 
spectrum of FPU lattice with periodic boundary conditions contains only two types of 
resonances. The first one is the 1:1 resonance between conjugate modes (k, N-k) and the 
second one is 2:1 resonance between vr-mode and the modes with indexes N/3 and 2N/3. 
Therefore, the right hand sides of equations [28] contain the terms similar to XmXnXj with 
indexes (m, n) 7^ (fc, N — k) and j=k or N-k. After integrating these equations over the 
fast time the requirement of the absence of any secular terms leads to that only the terms 
mentioned remain in the right hand sides. If we take an interest in a mutual dynamics of 
conjugate modes, one should choose the initial conditions for equations (1251) as follows: 

X,,i(r = 0) = 0, j^k,N-k. 

In such a case the right hand part items in equations ( !28l) are equal to zero and the 
resulting equations are uncoupled with any other modes: 

. , . = n 

(29) 



+ (/tilXfc.iP + /«2|XiV-fc,iP)Xfc,i + f^3XN-k,iX*k,i = 



+ (Ki|x7V-fc,lP + '^2|Xfc,lP)XJV-fe,l + l^3XllX*N-k,l 

_ 2,(3uk„ _ 3puJk„ _ ^pujk 

1^1 — — - — ^k,k,k,k, 1^2 — ^k,k,N-k,N-k, 1^3 — — 7, — ^k,k,N-k,N-k- 

o 4 o 

One can easily check that the coefficients Ck,k,k,k and Ck,k,N-k,N-k have the following 
values: 



4, k = N/A ' ^k,KN-k,N-k - j k = N/A - 

Such values of parameters show that the equations have the universal structure for 
any pair of conjugate modes. The only exceptions are the modes with indexes (N/4, 



3N/4] , which form the one-solution manifolds in the exact problem [ iPoggi and Ruffo. 



(119971 )]. Moreover, the following relation between coefficients in equations (1251) is valid: 

1^1 - l<'2 - 1^3 = 0, (30) 
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which is significant as we will see later. One can assert that this relation is valid for 
a/3-FPU lattices too. 

The Hamilton function corresponding to equations fl29|) : 



Hx = f (IXfel^ + IXiV-fcH + /«2|XfePIXJV-feP+ ^3-^^ 

^iXk XN-k + XkX*N-k) 

It is easily to show that equations (|29l) besides the integrals energy and momentum 
have an additional integral of motion 

Xk = \Xk\^ + \XN-k\'. (32) 

Here we discard the second term in the index of modes. The value is analogous 
to the quantum representation and it is often named as "occupation number" . It charac- 
terizes a full level of excitation of the conjugate modes. Taking into account the relation 
( I32l) one can rewrite the equations ( |29l) as follows: 



+ f^lXkXk + K3{xN-kXl - XkX*N~k)XN~k = 
+ f^l^kXN-k + KsiXkXll-k - XN-kX*k)Xk = 



(33) 



The first items in the equations (!33|) describe the change of eigenvalue Uk caused by 
excitation of the conjugate modes. The second terms describe the "exchange" interaction 
between modes. It is clear that equations (l33l) admit the one-mode solutions: 





If the initial conditions include the excitation of both modes, one can show that the 
value iG^ = {xN-kXl ~ X*N-kXk) is the integral of motion too. The energy of the system 
can be written as the function of Xk and G^: 

Ek - yX, - -G^ . 



Therefore the equations ( 1331) can be linearized: 

^fe + i^iXkXk - iKsG^XN-k = 

O 

The particular solution of equations (IMl) is 

Xk = akexp[i{QT2 + f 



+ f^iXkXN-k + t>^,G-Xk = ■ ^ ^ 



XN^k = a-N-k exp[z((]r2 - f )] ^^^^ 



with 



f2 — KiXk + li^G, 



and 



N-k 



A = ±V2- 
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This solution corresponds to the elhptical in-phase or anti-phase mode. The general 
solution can be presented in the form: 



XN-k = VJQe'^'^^^^ sin^e-^^/2 • 
The equations (IMl) in the terms of variables 6 and A have the following form: 

sm29[-^ + K3G,^ cos A] = 
sin2^[^ - 2K3G~ctg29smA] = ' 

At that, the parameter is connected with the variables (A, 9) by relation: 



(36) 



(37) 



sin 29 sin A = Gf. /X^ = const. 

The behavior of the system can be analyzed efficiently in the reduced phase plane 
(A, 9) (see Fig. [2]). 




Figure 7: Phase plane of equations ( !37l) in the terms (A, 9). The stationary point 
(A = ±7r/2, 9 = 7r/4) corresponds to elliptical mode (G^ = X^). The rectangular 
trajectories are the LPTs and they correspond to = 0. Both presented quadrants are 
equivalent ones because of periodicity conditions. 

There are two characteristic phase trajectories in Figd First of them is the central 
point (G^/Xfc = 1, A = 71/2, 9 = tt/4:) corresponding to elliptical normal modes fl35l) . The 
trajectories close to this point correspond to weak intermode energy exchange. Second 
trajectory (G^/X^. = 0) is limiting phase trajectory (LPT) describing complete energy 
transfer. Let us note that the rate of energy exchange rises with increasing Gl- and is 
minimal if G^/X^. = 0. So, both intensity and rate of energy transfer are controlled by 
the parameter G^/X^. The most intensive energy exchange is attained for the trajectories 
close to LPT (for LPT itself the rate is equal to ze ro that is st r ongly different from the 
case of two linearly coupled nonlinear oscillators [Manevit^h! ( 2007 )]. where complete 
interparticle energy exchange occurs with finite rate). 

Certain results of computer simulation of high frequency part of normal modes in 
the /3-FPU chain with 6 particles are shown in figure [HI 

Up to this time we have considered the dynamics of conjugate modes with initial 
phase shift 7 = 7r/4 (see eq. ffTI?]) ). One can easy to see that the harmonic part of Hamilton 
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Figure 8: Effect of energy exchange between conjugate modes (5,7) in the /3-FPU with 
12 particles, (a) The initial excitation corresponds to /X^ = 0. The absence of energy 
exchange between modes (left graph) leads to immobility of energy along the chain (right 
graph), (b) The initial excitation corresponds to /Xk ~ 0.1. The intensive energy 
exchange between high-frequency conjugate modes (left graph) leads to excitation mobility 
through the chain (right graph). 



function (fTSi) does not depend on the choice of 7 because any superposition of linear 
normal modes is a normal mode too. In the case of NNMs this statement is not more 
valid. Therefore, it is very interesting to estimate the effect of change of phase shift on the 
dynamics of NNMs. Let us consider the conjugate modes with indexes (N/4, 3N/4), th at 



form the one-mode solutions in the exact nonlinear problem Poggi and Ruffo. I (119971 )]. 
Therefore, the equations of nonlinear dynamics fl29l) do not contain any exchange terms, 
if the phase shift is equal to 7r/4: 

i-g:^ + i^i\Xn/4\ Xjv/4 = 

+ '^l|X3iV/4pX3Af/4 = 

For the /5-FPU lattices the phase shift 7 = which does not coincide with 7r/4, leads 
to the equations those are like the equations (IMIl . If 7 = is our choice for phase shift it 
is turns out that the relation 

Ki - K,2 + K,s = (38) 

is valid instead of eq.( l35i) . It leads to that the value 

is the integral of motion instead of the parameter G^ . In such a case the linearized 
form of equations ([33] )is: 

i-^^ + '^i^N/iXN/i + i^zG'^i^XzN/i = 
The stationary solutions have the following form: 
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X7V/4 = V^exp[z(nr2 + I)] 

X3N/. = ./^expmr, - f )] (40) 

with 

A = 0, G^i^ = Xn/a 

or 

A = 7r, G^y^ = — XAr/4. 

The main difference of these solutions from (135!) is that the displacement of stationary 
points from the positions (vr/4, 7r/2) and (vr/4, — vr/2) into (7r/4, 0) and (vr/4, vr) in the 
reduced phase plane (A, 6). The presence of exchange terms in the equations leads to 
that the trajectories, starting near the LPT describe the full energy exchange between 
modes. The comparison of dynamics of normal modes with different 7 is shown in the 
figS 
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Figure 9: The comparison of energy exchange between normal modes with different phase 
shift 7 in the /5-FPU chains with 4 particles, (a) The energies of one- and two-mode 
solutions at 7 = 7r/4 . There is no inter- mode energy exchange, (b) If 7 = the intensive 
energy exchange occurs at the same initial conditions. 

Equations (139|) in the terms of variables (A, 6) are written as 

sin2^[^ + «:3G+/,sinA] = 

sin 2e[j^ + 2K3G^/^ctg2e cos A] = ^ ^ 

The last equations determine a relation between angle variables in the form: 

sin 26^ cos A = = const. 

In this case the central point G^^jX^/i =1 corresponds to "supernormal mode" 
which is the combination of both resonant modes with identical weights. As for LPT 
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= 0), it describes again a complete intermode transfer with zero rate, so that 
actual slow energy transfer may be observed for the trajectories close to LPT. The data 
of computer simulation study confirm the analytical results (see fig. [TOl) . 




Figure 10: /5-FPU chain with 4 particles: (a) - (A, ^)-plane: the starting point is near 
the LPT; (b) the values X (black solid), (blue dashed) ,G~ (red dashed) versus time; 
(c,d) the time evolution of angles 6 and A (e,f) inter-mode energy exchange. 



The a/3-FPU chain is more complex object for analytical study. The structure of 
the phase plane depends on both parameter a and the number of particles in the chain. 
We use the a/?-FPU chain with 4 particles to clarify the main peculiarities of the system. 
The parameters of equations ( l33l) are written as follows: 

3/3 3/3 3/3 - 4a2 

«i = -r; ^2 = — ; ^3 



4 ' - 2 ' 4 
One can see that in such lattices the relations fl30f38l) are not more valid. It has two 
consequences: (1) both the parameter G~^ and the parameter G~ are not a constant of 
motion, and (2) the coefficient K3 can be equal to zero at a = ^3/3/4. All of that lead to 
more complicated description than discussed above. In such a case equations (15^ cannot 
be linearized, but they are written in the terms of angular variables as before: 

sin20[^ + i«:3Xsin2A]=O 

sin 20{f| -Xcos20[(/€i — K,2) ~ ^3 COS 2A]} = ' 
The first integral of equations (l42l) defines the relationship between variables (A, 6): 



sin 26*1/1^^ — — — cos2A| = C = const. 



K-3 

Here the phase plane has a structure depending on the relationship between the con- 
stants Ki, K,2, 1^3, which are controlled by the parameter a of the potential asymmetry. The 
case a=0 leads to relationship (|38l) with supernormal mode as stationary point (figlTOla). 
Any trajectory close to LPT corresponds to full energy exchange between modes xi ^"^^ 
X2 with large period. But occurrence of any asymmetry of potential function leads to 
separatrix creation. 

The singular (saddle) point for small values < K3 < ki — K2 (small a values) 
corresponds to unstable elliptical mode (A = 7r/2, = 7r/4). As a result we have got a 
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Figure 11: Transformation of the phase plane of equations fH2|) with a variation of poten- 
tial asymmetry parameter a (see text). 

phase plane shown in fig. (ITTl a). The separatrix crossing the point (A = 7r/2, Q = 7r/4) 
separates the domain of closed trajectories with partial energy exchange from the transit- 
time trajectories for which an energy exchange is practically absent. The domain of closed 
trajectories decreases while the parameter k,-^ tends to zero at a = 0.866 (fig. [TTlb). 
Two branches of separatrix are reduced to straight lines and the phase plane is filled 
with transit-time trajectories, which are mutually parallel ones. The supernormal mode 
(A = 0, ^ = 7r/4) appears as the saddle point within the interval -(ki — K2) < K3 < for 
0.866 <a < 1.2248 (fig. [He). And there are two stationary points (A = Tr/2, 9 = tt/A) 
and (A = 0, ^ = 7r/4) when the parameter K3 < (ki - K2) {a > 1.2248) (fig. [TTld). 
In such a case, the trajectories rounding both supernormal and elliptic stationary points 
and close to the LPT, correspond to intensive energy exchange between conjugate normal 
modes. So, the dynamics of the particles drastically changes with a variation of parameter 
K3, which, in turn, is determinate by the asymmetry of potential. 

To check the main analytical results we have performed the computer simulations 
of dynamics of several FPU-chains with various number of particles (N=3,. . . ,8). Some 
results of these simulations in the partial case N=4 are shown in the fig. [T21 The computer 
simulation was continued since t=0 up to t=1500, but the parameter of asymmetry was 
changed at t=300, 600, 900, and 1200. The starting points are inside the separatrix at 
a = 0.1 and a = 1.0, but it is in the transit-time area for a = 0.5. The last means that 
the energy exchange between modes is forbidden. Because separatrix is absent at « = 1.5 
the complete energy exchange corresponds any trajectories starting near the LPT. 

An example of energy exchange between conjugate high-frequency modes, which gives 
rise to the energy transfer along the chain is shown in the figures [T3]fT^ 
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Figure 12: Intermode energy exchange in the a^-FPU chain with 4 particles with various 
parameter a. The left column contains the phase trajectories in the (A, (9)-plane and the 
right column show the energy exchange between conjugate modes. 
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Figure 13: FPU chain with 6 particles: (a) Phase trajectory and angles A and Q] (b) 
energy of conjugated modes. 
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Figure 14: FPU chain with 6 particles: the energy transfer along the chain corresponds 
to the data shown in the figJTSl En - energy of the particles. 

Spatial energy localization in the discrete FPU chains. 

The dynamics of the normal modes in the essentially discrete system considered above 
has been studied in the assumption that the gap between vr-mode and nearby conjugate 
ones is large enough. But, if the number of particles grows the new opportunity arises. 
Namely, the resonant interaction between the modes having close frequencies (in linear 
approximation limit) becomes possible. Really, the frequency of 7r-mode differs from the 
frequency of nearby ones by the value of order l/N"^: 



TT /v TT 1 TT 

UN:, = 2sin(— (— ± 1)) = 2cos(— ) ^ujn(1- -( — ?). (43) 

Because the corresponding integral manifolds, two qualitatively different scenarios 
can be single out. For small enough energies a simple superposition of the normal vi- 
brations occurs that is similar to linear approximation. However, after attaining certain 
threshold the resonant interaction mentioned above occurs and it leads to some local- 
ization of excitation which can be move along the chain. This change in the behavior 
of the system can be interpreted as a manifestation of instability of invariant manifold 
in space of normal mode. The instability of zone-boundary vr-mode, well known in the 
FPU-chain s, usually corresponds wit h onset of weak chaos and the generation of chaotic 



breathers Lichtenberg et al. I (120081 )]. However, the results of the previous consideration 
show that the equation of dynamics of vr-mode can be effectively linearized. So nothing 
instability of vr-mode has been observed. Even if an initial excitation is distributed be- 
tween several modes (the vr-mode and nearby ones) the resulting dynamics of the chain 
is well described by simple superposition of these modes. The computer simulation data 
are in good accordance with the analytical model. 

While the energy of excitation and/or number of particles grows the dynamics of the 
chain becomes different than that for the quasi-harmonic dynamics defined by equations 
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(IT5]) . While the parameter e decreases, the eigenvalues of the vr-mode and the nearby 
ones become more and more close until the resonant conditions occur. 

In such a case the equations fl33|) are not valid because this resonance does not take 
into account. Now we choose the top frequency mode Xn/2 as the basic one. So, because 
the difference between Xn/2 and Xn/2-i has the order additional linear items appear in 
the equations for Xn/2-i and Xn/2+i- Moreover, we have to take into account the resonant 
interactions in the nonlinear part of equations. Missing out the details of calculations we 
get the final equations in the form: 



+ ^i[\XN/2?XN/2 + 2(|X7V/2-lP + 1x^/2+1^X^/2 + 

(XAr/2-1 + Xn/2+i)'X*N/2\ = 

^'-^ - TX^/2-1 + fmXN/2? + 3|XiV/2-lP + 

'^\XN/2+i\^)XN/2^1 + (2XAr/2 + XAr/2+l)XjV/2-l] 
'^Xjv/2+i vr^.^ , I 3/3r//iK^ , |2 I qU, , 12 



(44) 



I 



dT2 



'yXN/2+1 + f[{MXN/2? + 3|X7V/2+lP + 



'^\XN/2^i\^)XN/2+1 + (2X^/2 + Xlf/2~l)X*N/2+l] = 

First of all one should point out the presence of "exchange" items (x|XAr/2' 3 ~ 
N/2 + 1, N/2 — 1) in the equations (jBl) which are analogous ones to those in the equations 
for conjugate modes fl33l) . It is easy to check that equations (jH]) possess an additional 
integral of motion - total "occupation number" 



X=\Xo\ +IXi| +IX2I = const. 

One can see that the equations (jSj) admit both the one-mode and the two-mode 
solutions. If XAf/2 = is the choice we get the equations that are analogous of equations 
f l33|) . because the presence of the linear items does not introduce any significant distinc- 
tions. The computer simulation data confirm the existence of dynamical regime with only 
excited conjugate modes even in the long chains. On the other side, it is easy to show 
that a combination of modes with different eigenvalues (i.e. a combination of vr-mode and 
one of conjugate mode) leads to the localization of energy in a one half of the chain (see 
the fig. [ISD. 



Figure 15: The energy distribution in the /3-FPU system of 12 particles corresponding to 
the excitation of vr-mode in combination with one of conjugate mode. 

And at last, the energy put only in the mode with specified "parity" (e.g. Xn/2-i) 
is not transferred into mode with other "parity" (one can show that this statement is 
destructed by instability of conjugate modes). Let us consider the equations (jH]) with 
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initial conditions {xn/2+i = 0, Xn/2 0) Xn/2-i 0) ^or the determinacy. In such a case 
to analyze the dynamics of the system it is convenient to introduce new variables (po, ipi. 



XN,2 = ^M + ^^) (45) 

XN/2-1 = ^(<^0 - <^l) 

Because the "occupation number" X = |y)oP + Iv^iP is the integral of motion as 
before, the variables ipo and ipi can be expressed by the angle variables: 

The Hamilton function, written in the terms [Q^ A = (5i — (5o), has the form: 

R{B, A) = § [51/?X - 167r2 + 2(87r2 - 3/3X) cos A sin 20 

-3/3X(8 - cos2 A) sin^ 2^] ' ^ 

The respective equations of motion 



g + 3^ (Stt^ - + cos A sin 2Q) sin A = 

sin 2^11 - ^ cos 2^[(87r2 - 3/?X) cos A - 3/3X(8 - cos^ A) sin 2Q\ = (48) 
A = 6i- 6o 

have the stationary points: 

(a) A = 0, e=l 

(6) A = 7r, = f . (49) 

(c) A = 0, ^1 = larcsin(^^^Jj^), 9^ = ^ - 9^ 

Two first stationary points ((a) and (b)) correspond to equal values of |v9o| and \(pi\ 
, and it leads to the pure vr-mode (a) or pure XAf/2-i-niode (b) in the accordance with 
definition (HSil . These stationary states exist at every values of the amplitude X. The 
limiting phase trajectory (LPT) rounding the stationary point (a) corresponds to the 
transfer from the state where Xn/2 = Xn/2~i to the state where Xn/2 = —Xn/2-i, that 
leads to energy transition from one half of the chain to another one. That is the nonlinear 
analogue of beating phenomena in the linear system. In the small-amplitude limit it is 
the only mechanism of the energy redistribution in the chain, when the energy exchange 
between conjugate modes is absent. 

The stationary points (c) originate from the point (a) when the amplitude X reaches 
the value 7r^/3/5. At this instant the stationary point (a) becomes the saddle one, conse- 
quently the zone-boundary vr-mode loses its stability. The instability of vr-mode appears 
as a small admixture of nearby conjugate mode that leads to a weak modulation of am- 
plitude of original mode. If the amplitude X reaches large enough values, the modulation 
of vr-mode turns out be large, that leads to generation of short life-time breathers. The 
example of such instability in the chain with 50 particles is shown in the fig. dni The 
energy of the chain corresponding to the threshold value X is equal to vr^/3/?X (only own 
part of the mode en ergy is taken into accoun t). This result is in the good accordance 



with the estimation Lichtenberg et al. I (120081 )] in the limit of the large chain {N ^ 1 



While the amplitude X grows, the separatrix rounding the stationary points (c) 
enlarges but the transition along the LPT is possible up to the instant when the separatrix 
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Figure 16: Instability of 7r-mode in the /?-FPU chain with 50 particles: (a) 3D surface of 
energy of particles; (b) energy of normal modes. E - energy of particle, n - number of 
particle, t - time, Sk - energy of normal modes; the initial occupation number Xq = 6.25. 

merges with the LPT. Since this instant the trajectories starting from the point with 
6 < 7t/4 and = cannot reach any point with 6 > n/A. It means that the energy of 
initial excitations turns out to be confined at the one part of the chain. The threshold 
of such localization can be calculated from the requirement of the equality of energy 
of the states in the LPT and the energy of unstable stationary point (a). The result 
is Xc = 167r^/27/? that corresponds to the energy threshold Eioc = 167r^/27/?A^. It is 
convenient to observe the evolution of phase trajectories accompanied the variation of 
amplitude X on the phase space in the terms of angle variables A and 6. Four different 
structures of phase plane are shown in the fig. [T71 The first graph (ITTl a) shows the phase 
plane before the first bifurcation, i.e. at the small values of X. After the first bifurcation 
the small separatrix crossing the unstable point (0, 7r/4) originates at X = it'^/3i3 (fig. 
[T71b). The second bifurcation is shown in the fig. [T7] (c). It is well seen that the LPT 
merges with the separatrix. The subsequent rise of amplitude X leads to the creation of 
transit-time trajectory region (fig. [T7] (d)). 

The LPT starting the point (0, 0) is the boundary trajectory divided the area of 
closed motion around the stationary point (0, ^i) from the area of transit-time trajectories. 
The new separatrix crossing the point (0, 7r/4) describes the transition {xn/2 XN/2) and 
rounds the stationary point Xn/2-i- One should note that the LPT describes the periodical 
process which is accompanied by some deformation of energy distribution profile (the 
"breathing" mode of localized excitations - see figJT3lb). The period of this "breathing" 
can be estimated as 



where the last integral is taken along the LPT. The equation of LPT is defined by 
proposition 




(50) 



H{e,A) = H{0,0) 



X 



(51/?X- 167r2). 



64 



This relationship leads to 
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Figure 17: Phase plane transformation when the amphtude X grows. {ei)PX — n^/S 
Phase plane before the first bifurcation: there are only two singular points like center; (b) 
j3X = 7r^/2 after the first bifurcation the separatrix crossing the stationary point (0, 7r/4) 
describes the instability of 7r-mode; (c) the LPT merges with the separatrix at the critical 
value of occupation number Xc = 167r^/27/9, the transitions from any state below 9 = ir/A 
near ^ = to the state with 9 > n/A are forbidden and vice versa; (d) PX — 27r^/3 - the 
area of transit-time trajectories is observed at the any occupation number X > Xc and 
LPT is limited by the amphtude, corresponding to the angle 9 < 7r/4. 
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sm2e = 2k-, 



K 



"S-cos^ A 
8 7r^-3/3 X 
3/3X 



(51) 



The variation of period as a function of parameter k is shown in the figure [T8l 




8;r' -3/3X 

3;SX 



Figure 18: Period of motion along the LPT versus parameter k. The value k =3.5 
corresponds to the threshold of localization = 167r^/27/9. The smaller values of k 
respect to the localized regime with X > Xc- 

In the fig. [19] the energy distribution along the chain with 12 particles is shown both 
before threshold of localization and after that. 




(a) " (b) 

Figure 19: The energy distribution along the chain with 12 particles at X = 147r^/27/3 < 
Xc (a) and X = 27r^/3/3 > Xc (b) at the different times t. The "breathing" mode of the 
energy profile in the right side corresponds to motion along the LPT (see text). 



So one can see that energy localization in the one half of chain needs a simple super- 
position of two modes with the total occupation number X exceeding the critical value 
Xc- The value of Xc does not depend on the length of the chain but the total energy E is 
directly proportional eX = X/N that leads to the infinitesimal threshold of localization 
in the case of infinite chains. Really this process is an outset of the creation of localized 
excitations (such as chaotic breathers) in the system. The computer simulation data show 
the energy localization in a more narrow area if the chain is the long enough {N ~ 20 and 
more). Such a contraction of localization area needs a participation of modes different 
from the considered ones. The presence of these additional modes is reflected in the phase 
portrait of the system (flg. [201) where the computed phase trajectories are more compli- 
cated than the analytical predicted ones (flg. [T7|l . The flgure 15 shows the long time 
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localization in the system with 8 particles. One can see that the motion of the position of 
energy maximum correlates with energy exchange between conjugate modes Xn/2-i and 
Xn/2+1- The reason of that is the difference between asymptotic equations fj44|) and the 
exact dynamical behavior of the original FPU chain. These differences grow when the 
amplitude of excitations increase that manifests as an instability with respect to energy 
exchange between the conjugate modes considered. 




Figure 20: The LPT on the (6', A)-plane for the chain with 12 particles: (a) -before the 
threshold of localization and (b) - after that. These trajectories correspond to the energy 
distributions showed in the fig. [T9l 




Figure 21: The computer simulation data for the chain consisting 8 particles: (a) The 
surface of energy of particles, (b) the energy of normal modes. The initial excitation 
corresponds to X = 8.0 . 



Conclusion. 

The dynamics of essentially discrete nonlinear systems with periodic boundary conditions 
demonstrates two typical regimes of behavior. The first of them is due to the discrete- 
ness of the spectrum of eigenvalues of the system and corresponds to the quasiperiodic 
motion along the invariant tori in the Fourier-space of the system. In such a case the 
dynamics of resonant pairs of conjugate modes in small FPU systems can be described 
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from unique viewpoint using complex representation of dependent variables and multiple 
scale techniques. As this takes place, the complete general classification of possible types 
of dynamical behavior is attained. The most simple scenario is reahzed for the case of 
symmetric nonlinearity (/3-FPU system) when only two characteristic phase trajectories 
corresponding to supernormal (elliptical) mode and LPT, can be discriminated. They 
describe the vibrations with phase shift A = (A = 7r/2) between two resonant modes 
and complete energy transfer between them, respectively. Both intensity and rate of in- 
termode energy exchange are controlled (depending on the resonant pair) by parameters 
Gf /Xi or Gl /Xi (despite the difference between Gf /Xi and G^ /Xi, corresponding sce- 
nario are similar except of the presence of supernormal mode in the first case and elliptical 
mode in the second one). Rise of this parameters leads to decreasing of exchange intensity 
and increasing of its rate. It is very important that the evolution of initial excitation is 
described efficiently by the set of quasihnear differential equations parametrized by the 
"occupation number" of the respective normal modes X. 

Introduction of asymmetry leads to complication of the dynamics of the conjugate 
modes with indexes iV/4 and 3iV/4 with several qualitative transitions. The evolution of 
dynamical behavior with increasing of asymmetry is inverse (in certain sense) to that for 
linearly coupled nonlinear oscillators. For strong asymmetry two stable stationary states 
corresponding to supernormal and elliptical modes exist. Two branches of LPT describe 
complete energy exchange. When parameter of asymmetry decreases the instability of 
supernormal mode occurs with formation of singular trajectory (separatrix) and breaking 
off the complete energy transfer. Further evolution with decreasing of the parameter of 
asymmetry leads to full termination of energy exchange between resonant modes, restor- 
ing the stability of supernormal mode with simultaneous instability of elliptical one and 
possibility of partial intermode energy exchange. If asymmetry disappears the separatrix 
transform to LPT with complete energy transfer. 

The second dynamic regime is associated with the interaction between non-conjugate 
high-frequency modes, i.e. modes with different eigenvalues in the linear limiting case. 
However, the small splitting of the frequencies of zone-boundary 7r-mode and modes 
nearby of it provides almost resonant conditions between them. This circumstance can be 
the consequence of both elongation of chain length and the rise of the amplitude of oscil- 
lations. It is very interesting that the natural small parameter here is the inverse length 
of the chain 1/N. The quasi-resonant interaction is a reason of instability of 7r-mode when 
the energy of the chain reaches the first threshold value E — n'^/SPN. This energy cor- 
responds to the occupation number X — 7r^/3/3, which does not depend on the length of 
the chain. The subsequent rise of occupation number up to value X^ = 167r^/27/9 leads 
to the confinement of the energy in the local area of the chain. The reason is that the 
separatrix crossing the unstable stationary point corresponding the 7r-mode increases and 
merges with LPT. This circumstance leads to the breaking of trajectories corresponding 
to the energy transfer from one half of the chain to other one. The localized excitation 
corresponds to new LPT demarcates the areas of closed trajectories and transit-time one. 
The energy profile of localized excitation undergoes some disturbances while the imaging 
point of the system moves along the LPT. It is a cause of "breathing" mode of the chaotic 
breathers a starting point of which is the process of localization discussed above. Moreover 
increasing of the amplitude of oscillations in the nonlinear system leads to the frequency 
shift that makes the interaction of high-frequency modes with more low frequency ones. 
The computer simulation data show some contraction of area of energy localization ac- 
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companied by the excitation of low-frequency modes. So the chaotic breathers result to 
the interaction of many modes with different eigenvalues regardless of the fact that the 
contribution of low-frequency modes can be small enough. The last remark concerns with 
the role of boundary conditions. In spite of the fact that the periodicity of the system 
studied is very essential proposition the systems with other boundary conditions can be 
considered in the framework of this approach. The decreasing of frequency splitting of 
the top modes with different eigenvalues occurs while both the length of the chain and the 
amphtude of oscillations increase. Therefore the efficient resonant interaction can result 
in the localization and confinement of energy in the bounded area of the chain. 
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